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A theoretical and numerical analysis of the transition from chaotic to 
nonchaotic behavior in an ensemble of particles with different initial 
conditions which move according to Newton's equations in a bounding 
potential and are driven by an identical sequence of random forces (see 
S. Fahy and D. R. Hamann, Phys. Rev. Lett. 69, 761 (1992)) is 
presented. The threshold values of the parameters for transition from 
chaotic to nonchaotic behavior are defined on the basis of the map for 
distances between the particles and differences of velocity. Numerical 
analysis is fulfilled for one-dimensional Duffing V(x) = x A — x 2 and 
V(x) = x 4 potentials. 
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Recently an interesting transition from chaotic to nonchaotic behavior 
in randomly driven systems has been discovered [1]. When an ensemble of 
bounded in a fixed external potential particles with different initial conditions 
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are driven by an identical sequence of random forces, the ensemble of trajec- 
tories may become identical at long times. Fahy and Hamann [1] considered 
a particle of mass m moving according to Newton's equations in a potential 
V(x), except that at regular time intervals r the particle is stopped and its 
velocity is reset to random value chosen from a Maxwell distribution with 
temperature T. It should be stressed that for every particle of the ensemble 
it was given an identical, randomly chosen velocity at the start of each step of 
time length r. This motion is in many respects similar to Brownian motion 
of the particles at a temperature T and is frequently used in Monte Carlo 
simulations [2]. However, if the time interval r between stops is lower than 
a threshold value r c , the final trajectories of the particles are independent 
on the initial conditions; all trajectories become point by point identical in 
time. Although the trajectory is highly erratic and random, the system is 
not chaotic. 

The similar effect may also be observed in a more general and realistic 
(from the physical point of view) case, i.e., when mixing at time intervals r 
some part a of the old velocity v old with random velocity v mn to get a new 
starting velocity v™ = av oW + v mn [3] . Here a threshold value r c depends 
on a. 

Note, that the observed effect resembles a phase transition but does 
not depend crucially on the dimension of the space in which the particles 
move. Besides its intrinsic interest, this phenomenon has some importance 
for Monte Carlo simulations [1] and, probably, can have influence on some 
processes, e.g., clustering of particles. 

In the paper [1] except the discovery in the computer simulations of the 
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transition from chaotic to nonchaotic behavior, the analytical estimation of 
the average rate of contraction of the distance between initially close together 
identically driven particles (for a one-dimensional potential, very short time 
intervals r between stops and a = 0) is given. Note to the paper [4] in 
which the similar problem was analysed and it was found that trajectories 
of a noisy logistic map for various initial conditions may become point by 
point identical in time if the noise is sufficiently strong and independent of 
the initial condition. Moreover, using Langevin equation, it was shown quite 
generally, that nonlinear systems coupled through identical strong noises may 
have the same trajectory in the long time limit. However, the derivation of 
dependence of the threshold noise on the parameters of the system [1] is 
absent. 

It is the purpose of this paper to present a theoretical and numerical 
analysis of the transition from chaotic to nonchaotic behavior in randomly 
driven one-dimensional systems, i.e. the dependencies of the threshold value 
r c on the parameters of the problem, including a, and convergence rates of 
the trajectories. Generalization of the method for two- and three-dimensional 
systems is rather straightforward (see below). It should be noted that in 
the numerical calculations of the mean-square distance between particles we 
observe the intermittency route to chaos when r — > r c : mean-square distance 
as a function of time alternates randomly between almost clustered phases 
and irregular bursts of the dispersed particles. Moreover, we find a transition 
from chaotic to nonchaotic behavior not only for the regular time intervals 
t but for the randomly changing intervals between stops, i.e. for random Tj 
(where % is the number of the step) as well. Our theoretical analysis is based 
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on the mapping form of the equations of motion for the distance between the 
particles and the difference of the velocity of the particles. 

Let us consider a particle of mass m moving in a one-dimensional potential 
V(x) which confines particles to a finite region. At a time intervals Tj the 
particle is partially stopped and its velocity is reset to a new starting velocity 
Vi = avf dj rv\ an . Between the stops the particle moves according to Newton's 
equations 

dx d 2 x 1 dV . . 

dt dt 2 m dx 

When two particles initially at points x$ and x' are started with velocities Vq 
and v' and are driven by an identical sequence of random velocities v \ an at the 
same time intervals Tj, coordinates and velocities of them may accidentally 
draw closer to one another. The convergence of the two trajectories to the 
single final trajectory will depend on the evolution with a time of the small 
variancies of the distance Axi = x[ — Xi and velocity Avi = v\ — Vi. Moreover, 
we investigate a transition from chaotic to nonchaotic behavior. Generally, 
such a transition may be detected from analysis of behavior of the neighboring 
trajectories and it is described by the Lyapunov characteristic exponents and 
KS metric entropy of the flow of trajectories in a given region of phase space 
[5-9]. Further we will follow this method. Comparisons of the results of such 
analysis with investigations by the direct numerical simulations confirm the 
suitability of our theoretical approach. 

From formal solutions t) and v — v(xi,Vi,t) of equations (1) 

with initial conditions x = Xi and v — at t — it follows an equation for 
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Ax(t) and Av(t) at a time moment t: 

( Ax{t)\ 



(£w) =T(a:w) 



(2) 



where the matrix T is of the form 



T xx (%T XV 



/ dx(xj,Vi,t) dx(xi,Vj,t) \ 

dxi dvi 

dv(xi,Vi,t) dv(xi,Vi,t) 

V Fir,; r)ih ) 



(3) 



dxi dvi 

Note, that the similar method of investigation is used in the theory of transi- 
tion to chaos in classical systems [5,6]. However, the motion in the form (2) 
and (3) is represented as the non-area-preserving tangent map, while clas- 
sical dynamics of the conservative systems may be represented by the area 
-preserving maps. 

According to equations (1) and (3) matrix elements T xx and T xv satisfy 
the equation 

^2 T i J2 V 

T x (4) 



dt 2 



m 



dx 2 



X=x(xi,Vi,t) 

while T vx = T xx , T vv = T xv and the initial conditions at t — are: 



T xx {xi ) fj, 0) T vv 1, T xv T vx 

1 d 2 V 



(5) 



T xx (^Xi. ) Vi, 0) T vv 0, T xv 1, T v . 



m dx 2 

Therefore, the dynamics of the distance between the particles Ax and the 
difference of the velocity Av may be represented by the non-area-preserving 
mapping form of the equations of motion 



/ Ax i+1 \ , ( Axi 



(6) 



In general, the intervals between the resets of the velocity Tj may be depend- 
ing on the number of step i. 

Further analysis of the model may be based on the general theory of the 
dynamics of classical systems represented as maps [5-9]. Thus, for a = the 
Lyapunov exponent is defined as 

1 N 1 

A = Jim — V — In T xx (xi, v t , n) (7) 

and may be easily evaluated numerically 

For a — 1 the map (6) is area-preserving and det T(l; x iy Vi, 7$) = 1, while 
in general det T = a, TrT = T xx + aT vv and the eigenvalues /ii j2 of the T 
matrix are given by equation 

fi 2 - yuTrT + det T = (8) 

which yields 

1 



^1,2 = 2 



T xx + aT vv t \J (T xx + uT vv f - 4a] . (9) 
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So, the eigenvalues come in reciprocal pair, //1//2 = ot. For (T xx + otT vv y 
4a < the eigenvalues form a complex conjugate pair with = \/i 2 \ = 
otherwise the eigenvalues are real. 

Among the randomly driven systems only dynamics of the harmonic oscil- 
lator with regular time intervals r for the reset of velocity may be described 
exactly. In this case the T matrix is independent on the starting coordinate 
Xi and velocity v iy i.e. 



cos out (a/uj)smujr 
— uj sin lot a cos lot 



(10) 



If a < 1 and cos cur 7^ =pl modules of both eigenvalues of the T matrix 
are < 1- So there is a contraction of the distance in the phase space 

between the particles in each step. The transition matrix for n steps T n = T n 
may be calculated by means of diagonalisation of the matrix (10). The 
eigenvalues of the T n matrix are /z™ 2 and the Lyapunov exponents are simple 
Ai )2 = 7 In |/xi, 2 | - 

Generally, the mapping T(a; Xj, i>j, Tj) in (6) is depending on the starting 
coordinates and Uj. Therefore, calculation of the mapping for n steps, T n = 
T(o;; Xj +n _i, i>j + „_i, Tj +n _i) • T(o;; Xj + „_2, fj +n „2j r i+n-2) ■ ■ ■ T(a; Xj, f j, Tj), 
and of the corresponding eigenvalues are complicated problems. Further we 
will evaluate the averaged quantities 




which are analogous of the averaged Lyapunov exponent (7), characterise the 
rate of the exponential increase of the separation of the two initially adjacent 
points and are related with the KS entropy of the system [5,6]. Comparisons 
of the threshold values r c from the direct numerical simulations with those 
from the criterion 

^largest = (12) 

indicate to the usefulness of the quantities (11) for analysis of transition from 
nonchaotic to chaotic behavior of the systems. 

In Fig. 1 we show typical results of the direct numerical calculations 
of the mean-square distance between pairs of identically driven particles in 
the Duffing potential V(x) = x 4 — x 2 vs number of steps for different time 
intervals r between stops and a — 0. The averaged characteristic exponents 
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as functions of time between stops, r, calculated according to Eq. (7) for 
Duffing and V(x) = x 4 potentials are shown in Fig. 2. Regions of r with 
the positive and negative characteristic exponents correspond respectively 
to the chaotic and nonchaotic behaviors. Such a conclusion coincides with 
the conclusions from the direct numerical investigation of the systems. In 
Fig. 3, we show dependencies of the threshold values for transition from 
chaotic to nonchaotic behavior of the time r c between resets of the velocity 
to v ™ ew = avf d + v \ an on the parameter a calculated according to criterion 
(11)-(12). 

For the relatively short time intervals one can expand the position and 
velocity of the particles in the powers of time t. Using equations of motion 
(1) and relations (3) we have expansions for the matrix elements T xx and T vv : 

Txx = I - J-f^t 2 - ^—t 3 + ■■■ 
2m dx 2 6m dxf 

T I 1 d 2 V t2 Vj d 3 V^ | 
2m dx 2 2>m dxf 

We can use expressions (13) for the analytical evaluation of the average 
rate of contraction of the distance between trajectories for small r. So, the 
result of expansion of eq. (7) according to eq. (13) in powers of r is 

\r = -— I^X\ T 2 + ... 
2m \ dx 2 I 

where the angle brackets denote the average with respect to x,. After the 
averagement according to the Boltzmann distribution with ksT = m = 1 
(instead of the average with respect to X{ from random trajectories) we have 
Ar ~ — 2.2r 2 and Ar ~ —2. Or 2 for the Duffing potential and V(x) = x A 
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potential, respectively. For r < 0.5 these evaluations of the rate coincide 
with the numerical calculations in Fig. 2. 

Moreover, expressions (13) may be used for the description of the systems 
with a — > 1 and for short time intervals r. Introducing the notation a = 
1 — 2s (e — > 0) one can have from equations (9) and (11) the approximate 
expressions 

fii,2 — 1 — £ =F ^T xx + T vv -2 + e 2 + ■■■ (14) 

<7i,2 — - \— £ T ((T ra + T OT - 2 + e 2 ) 1 ' 2 )} • (15) 

e = (1 - a)/2 < 1, r < r c . 

The angle brackets in eq. (15) denote the average with respect to Xi and Vi 
only if T xx + T vv — 2 + e 2 > 0, otherwise /ii j2 are complex and cr 1 ^ — —e/r. 
The threshold value r c of time interval between steps for transition from 
nonchaotic to chaotic behavior may be evaluated from the criterion (12). 

Using the average according to the Boltzmann and Maxwell distributions 
(instead of the average with respect to and v « from random trajectories) 
one can investigate the behavior of the system according to equations (12)- 
(15). This yields to the expressions of the threshold time intervals 

r c ~ 1.44(1 - a )/[l + 1.6(1 -a)}, (16) 

r c ~ (1 - a) 2/3 , 0.01 < 1 - a < 1 (17) 

for the Duffing and V(x) = x 4 potentials, respectively. Here and in the 
numerical simulations we set k^T = m = 1. It should be noted, that the 
"effective temperature" of the system in such simulations is T/(l — a 2 ). 
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In Fig. 3 the dependencies (16) and (17) are compared with those ob- 
tained from the numerical calculations according to equations (4), (9), (11) 
and (12). Note to the dependence of the law of the approach r c — > as a — > 1 
on the shape of potential. 

All such numerical analysis and comparisons with the theoretical results 
indicate to the usefulness and fitness of the present theory for the investi- 
gation of transition from chaotic to nonchaotic behavior in randomly driven 
systems. Generalization of the theory to more degrees of freedom is straight- 
forward: for two- or three-dimensional systems equations (l)-(4) and initial 
conditions (5) for the T matrix are essentially the same, only instead of one 
equation (4) we have the systems of two or three coupled equations and di- 
mension of the T matrix is 4 or 6, respectively. Further analysis is rather 
similar to that according to equations (6)-(12). The only additional compli- 
cation is the diagonalisation of the multidimensional T matrix. 

In summary, we have developed a theory for transition to nonchaotic be- 
havior in a Brownian-type motion. Theoretical analysis based on the map- 
ping form of equations of motion for the distance between the particles and 
the difference of the velocity allows to simplify the problem of investigation 
of transition to nonchaotic behavior and results to the expressions for the 
criteria of the nonchaotic motion. Theoretical results agree well with the di- 
rect numerical simulations and indicate to the possibilities of generalization 
of the model, e.g. to more degrees of freedom, for random values of the time 
intervals between the resets of the velocity and for systems driven by the 
random forces. 
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Fig. 1. The mean-square distance (Ax 2 ) between pairs of identically 
driven particles in the Duffing potential V(x) = x A — x 2 vs number of steps 
for different time intervals r between stops. 500 independent simulations 
with initial conditions from the Gaussian distribution of variance a 2 = 1 are 
averaged. The time r between stops was r = 1.75 for the first 300 steps and 
r = 1.3 for the next steps. 

Fig. 2. Averaged characteristic exponents (Lyapunov exponents multi- 
plied by r) vs the time r between stops according to equation (7) for (a) 
Duffing potential and (b) V(x) = x 4 potential. 

Fig. 3. Dependencies of the threshold value for transition to chaotic 
behavior of the time r c between resets of the velocity to v™ ew = avf d + v r i an 
on the parameter a according to criterion (11) - (12) and according to the 
analytical limiting expressions (16) and (17) for (a) Duffing potential and (b) 
V(x) = x 4 potential. 
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